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We report a search for the decays B° — ► DJ" Df , B° — > D* s ~ Df and B° —* D* s ~ D* s + in a 
sample of 232 million T(4S) decays to BB pairs collected with the BABAR detector at the PEP- 
II asymmetric-energy e + e~ storage ring. We find no significant signal and set upper bounds for 
the branching fractions: B{B° -> D~ D+) < 1.0 x 10~ 4 , B(B° -> D*~D+) < 1.3 x 10~ 4 and 
B{B° -> D* s -D*+) < 2.4 x 1(T 4 at 90% confidence level. 
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In the Standard Model (SM), B° -» D y s 
cays are dominated by the W-exchange mechanism bd — > 
cc as shown in Figure 1, while the corresponding loop 
diagram is highly suppressed. The decay rates of W- 
exchange or annihilation processes are usually argued to 
be negligibly small due to the suppression from helicity 
and/or form factors [1]; however this assumption has not 
been well tested experimentally. 
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FIG. 1: W-exchange decay diagram for B° — > Ds D. 

Recently, it has been pointed out that it is difficult to 
calculate these decay amplitudes using the factorization 
approach, and a perturbative QCD (pQCD) [2] model 
has been used to predict the branching fraction for these 
decays. Prediction of branching fractions from an alter- 
native model [3] gives an estimate of non-factorizable con- 
tributions coming from chiral loops (CL) and tree level 
amplitudes generated by soft gluon emission forming a 
gluon condensate (GC) and it differs from pQCD ap- 
proach by large amounts, as shown in Table I. 



TABLE I: Summary of theoretical predictions of the branch- 
ing fractions. 

B Decays Branching Fraction (xl0~ b ) 

pQCD [2] CL-GC [3] 



B° -> D~D+ 
B° -> D* S ~D+ 
■ Dt-D* s + 



3° 
B° 



7-8±?;g 
6.0±if 

8.5±? : g 



25.0 
33.0 
54.0 



It has been estimated that a CP asymmetry of the 
order of 10% could arise between B° — ► DjDf and its 
charge conjugate [4] . A measurement of the decay rates of 
B° -> Di* ] ~Di* )+ relative to those of B° -> £)(*)-£>«+ 
will provide an estimate of the W-exchange contribution 



to the latter decay, a crucial piece of information for ex- 
tracting the CKM angle 7 from B° -> £>(*)-£>(*)+ an d 
B° -> Z)(*)-L>i* )+ decays [6]. 

Using 211 fb _1 of data taken on the T(AS) reso- 
nance with the BABAR detector at the PEP-II asym- 
metric B factory, we report a search for B° — > DjDf, 
B° -> D*~D+ and B° -» D*~D*+ decays [7]. We 
use the D*~ decays into and D~ decays into <fin~ , 
K° S K-, and K*°K~. The 0, K s and K*° mesons are re- 
constructed in their decays to K + K~ , tt + tt~ and K + n~ , 
respectively. 

The BABAR detector is described in detail elsewhere [8] . 
Tracking of charged particles is provided by a five-layer 
silicon vertex tracker (SVT) and a 40-layer drift cham- 
ber (DCH). Discrimination between charged pions and 
kaons relies upon ionization energy loss (dE/dx) in the 
DCH and SVT, and upon Chcrcnkov photons detected 
in a ring-imaging detector (DIRC). An electromagnetic 
calorimeter (EMC), consisting of 6580 thallium-doped 
Csl crystals, is used to identify electrons and photons. 
These detector sub-systems are mounted inside a 1.5-T 
solenoidal superconducting magnet. Finally, the instru- 
mented flux return of the magnet allows us to discrimi- 
nate muons from other particles. We use the GEANT4 
Monte Carlo (MC) [9] program to simulate the response 
of the detector, taking into account the varying acceler- 
ator and detector conditions. 

Charged tracks used in the reconstruction of 4>, 
K* (892)° and D s meson candidates must have a dis- 
tance of closest approach to the interaction point of less 
than 1.5 cm in the transverse plane and less than 10 
cm along the beam axis. All kaon candidates must pass 
particle identification (PID) criteria, based on a neural- 
network algorithm which uses measurements of dE/dx 
in the DCH and the SVT, Chcrcnkov angles and the 
number of Cherenkov photons in the DIRC. No PID 
requirement is applied to the pion candidates. A (f> 
candidate is composed of two identified kaons of oppo- 
site charge that are consistent with originating from a 
common vertex. We accept (f> candidates with invari- 
ant mass 1.000 < m K+K - < 1.039 GeV. K° can- 
didates are composed of two oppositely-charged tracks 
coming from a common vertex with an invariant mass 
0.491 < m 7r+7r - < 0.505 GeV. K* (892)° candidates are 
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reconstructed from two oppositely-charged tracks, where 
one track is identified as a kaon, with an invariant mass 
0.842 < m K - n + < 0.942 GeV. 

We reconstruct D~ mesons from decays to 07r~ , 
K®K~, and K*(892)°K~ using tracks coming from a 
common vertex with a x 2 probability greater than 0.1%. 
The reconstructed mass of Dj candidates is required to 
be within 2.7 to 3.0 standard deviations of the nom- 
inal mass; a typical mass resolution of D s is about 
5.1 MeV. The selected D 3 candidates are then kinemat- 
ically fit with their masses constrained to the nominal 
value [10]. In the decays Dj 0tt~ {K*(892)°K-), 
the <j> (K* (892)°) mesons are polarized longitudinally. 
Therefore the cosine of the decay angle Oh between the 
direction of the K~ from <j) from K* (892)°) and the 
Dj direction in the <j> (K* (892)°) rest frame is expected 
to follow cos 2 9h distribution. Background events from 
random combinations are expected to be uniformly dis- 
tributed in cos9h- We place a decay mode-dependent 
requirement on the minimum value of | cos Oh \ > which 
varies from 0.3 to 0.5 and rejects 13 to 24% of the com- 
binatorial background. 

D*~ candidates are formed by combining D~ and 7 
candidates with a mass difference AM = M D »- — M D - 
in the range of 0.125 < AM < 0.160 GeV. The photon 
energy measured in the EMC is required to be more than 
100 MeV. 

B° meson candidates are reconstructed by combining 
either (i) two oppositely charged D s candidates, (ii) one 
D* candidate and an oppositely charged D s candidate or 
(iii) two oppositely charged D* candidates. Finally, two 
quantities are used to discriminate between _B°-meson 
signal and background: the beam-energy-substituted 
mass niEs = -J E* } 2 — (pjjj) 2 and the energy difference 
AE = E* B - El, where El is the beam energy in the 
center of mass (CM) frame, and {E* B ) is the CM mo- 
mentum (energy) of the _B°-meson candidate. For signal 
events rn.ES peaks at the J3°-meson mass with a typical 
resolution of 2.5 MeV, dominated by the uncertainty of 
the beam energy, and AE peaks near zero indicating that 
the B decay candidate has a total energy consistent with 
the beam energy in the CM frame. Depending on the 
particular B° decay mode, the measured resolution for 
AE is 6.5 - 13.3 MeV. 

Multiple candidates are found in 3% to 5% of the se- 
lected events in the three different B° decay modes. The 
best candidate in each event is selected based on the 
smallest x 2 combination, where 



E 



AM- AM 



a AM 



(1) 



( I 

and the sum is over D y s and D y s candidates partici- 
pating in a particular B decay. The mean values (rn D ± 
and AM) are the nominal values given in Ref. [10] and 



the errors (a m ± , a am) are measured in a data control 

D s 

sample of B v — > D~ D s decays. 

A small source of remaining background is e + e~ — > qq 
production, which is suppressed based on event topology. 
We restrict the angle (9t) between the thrust axis [11] of 
the B° meson candidate and the thrust axis of the rest 
of the particles in the event. In the CM frame, BB pairs 
are produced approximately at rest and form a nearly 
uniform distribution in |cos0t|- In contrast, hadrons in 
qq events are produced back-to-back in two jets, which 
results in a | cos 0t | distribution peaked at 1 . Based on 
the background level of each mode, we require the value 
of I cos 9t\ to be less than a mode-dependent upper limit, 
which ranges from 0.83 and 0.9. We require R2 < 0.4, 
where R2 is the ratio of the second Fox- Wolfram mo- 
ment to the zeroth moment [12], both determined us- 
ing charged tracks and unmatched neutral showers in the 
event. 

For different B° meson decays, a signal region is de- 
fined in a two dimensional scatter plane of toes and AE 
as shown in Table II. Optimization of the selection is 
performed separately for each of the three B° decays [13] 
by maximizing a figure of merit, S 2 /(S + B), where S 
is the number of signal events in the signal box as de- 
rived from the MC simulation and B is the number of 
background events estimated from simulations of generic 
.B-decays and qq continuum. We use the same selection 
criteria for different B° decay modes if the figure of merit 
differs by less than 10%. 

After the aforementioned selection, four possible back- 
ground sources are considered. First, the amount of com- 
binatorial background in the signal region is estimated 
from the grand sideband region: —0.25 < AE < 0.25 
GeV and 5.20 < m ES < 5.27 GeV. The second source 
of backgrounds arises from B meson decays such as 
B° -> D^-D ( * )+ and B~ D { * ] ~ . These back- 
ground events have the same toes distribution as the 
signal, but their reconstructed energy is higher than the 
beam energy. Third, the cross-feed background that may 
arise among the six combinations of D S D S modes and the 
three reconstructed B° decay mode was studied with a 
large sample of signal MC and the corresponding contri- 
butions were found to be small. Finally, rare B decays 
into the same final state particles, such as non-resonant 
B° —> D~ K°K + , have the same Toes and AE distribu- 
tions as the signal. This source of background is esti- 
mated to be negligible. 

Figure 2 shows the distributions of candidates for (i) 
B° -» D-D+ (ii) B° -» D*~D+ and (iii) B° -» D*~ D*+ 
decays in the AE versus toes plane after all selection 
criteria have been applied. We find 6, 4 and 3 candi- 
date events in the signal boxes that survived the selec- 
tion criteria for the B° -» Dj D+ , B° -> D*~D+ and 
B° — > D*~D* + decay processes, respectively. The com- 
binatorial background in the signal box (A^£™ b ), is esti- 
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FIG. 2: Distributions of events in the A_E versus tties plane, 
for (i) 5° -> DJD+ (ii) B° -> D* S ~D+ and (iii) 5° -> 
Dg~Dg + decays after all selection criteria are applied. The 
box in each plot is the signal region based on studies with MC 
simulation as described in the text. 



mated from the number of events in the grand sideband 
region of the data. We compute the average number of 
background (N^) within the region E-x < AE < E\ 
GeV and 5.20 < m E s < 5.27 GeV from a fit to the AE 
distribution of the data events in the grand sideband (de- 
scribed well by a first order polynomial function P(AE)) 

P(AE), where 



as 



AT avg 
ly bkg 



N GSB x jj P(AE)/ r_ Q r> 



r-0.25 



.El and Ei are the AE energy bounds of the signal box 
as shown in Table II and Nqsb is the total number of 
events in the grand sideband region. The vtles projection 



of these background events is modeled with the threshold 
function [14], 



dN 
dx 



r^/l-x 2 /E£ 2 exp[£(l 



(2) 



characterized by the shape parameter £, the endpoint pa- 
rameter Et fixed at 5.289 GeV and x = ities- Ni 



bkg 



m 



the signal box is then estimated from scaled by a 



MeV to 

com 
bkg 



factor: / 5 5 2 2 7 9 ™/ J 5 5 f We vary E* b by ±2 
include its effect in the systematic uncertainties in N££™ b 
The measured uncertainties due to the choice of thresh- 
old parameter £, endpoint parameter E£, and parameter 
of the polynomial fit are combined in quadrature with the 
Poisson fluctuation of the number of events in the grand 
sideband to obtain the total error on iV^™ h . This pro- 
cedure does not account for any potential backgrounds 
that are enhanced in the signal region. The simulation 
indicates that only a small component of the background 
from the B° -► d1* )_ £)(*)+ and B° -> D<* )- £)W a decay 
exhibits a peaking toes distribution. This component, 
^bkg k ' i s extracted from a binned likelihood fit to the 
toes distribution of simulated events using a combination 
of the threshold function and a Gaussian. The £ param- 
eter in the threshold function is fixed to the value we ob- 
tained from the fit to the data grand sideband. The mean 
and width of the Gaussian component is fixed to the fit 
values obtained from B° — > D~ D + decays in the data. 
Uncertainties in arising from the Di* ] ~ and 

£)(*)- j)(*)o branching fractions [10] are added to its sta- 
tistical error obtained from the fit. V^?™ 6 and N^ k are 



bkg 



added to obtain the total estimated background, Nbkg, 
as quoted in Table II. 



TABLE II: The number of signal candidates (N ca nd), total 
estimated background (Nbkg), efficiency from MC simulation 
times the branching fraction (ej x B), and 90% C.L. upper 
limit for B° -> D~Dt, B° -» D*~ D+ and B° -> D*~D*+ 
decay modes. 





B° D~D+ 


B° -► D* s ~ D+ 


B° -> D* s - D*+ 


AE(MeV) 


-18.0-18.0 


-25.0-20.0 


-46.0-30.0 


m ES (GeV) 


5.27-5.29 


5.27-5.29 


5.27-5.29 




3.51 xlO~ 4 


1.47 xl0~ 4 


0.85 xl0~ 4 


N cand 


6 


4 


3 


Nbkg 


3.3 ±1.0 


3.9 ± 1.2 


2.3 ±0.9 


U.L. 


< 1.0 x 10" 4 


< 1.3 x 10~ 4 


< 2.4 x 10" 4 



We consider the following sources of systematic un- 
certainty for the signal efficiencies. The particle recon- 
struction and identification efficiencies are obtained from 
simulation, and cross-checked and corrected using large 
data control samples. This results in systematic uncer- 
tainties of (1) 0.8% per charged track; (2) 2.5% per recon- 
structed K® candidate; (3) 2.5% per identified charged 
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kaon and (4) 1.8% per reconstructed photon. The uncer- 
tainty on the number of BB events is estimated to be 
1.1%. Depending on the B submodes, the error from 
the MC statistics is 2% to 4.5%. The systematic er- 
rors are dominated by the 13.3% relative uncertainty 
on B(D- -> 4>-k~) [15], and 15.8% and 9.8% errors in 
B(D- -> K S K~) and B(D~ -> K* a K~) relative to 
B(D~ — > 07r~), respectively [10]. The uncertainty in 
modeling the simulation of the AE, | cos 9t | , | cos 9r \ 
distributions is evaluated using a ratio of the signal yield 
from B° — ► D~Di*^ + data control sample and generic 
BB MC. Each selection requirement is varied and the 
resulting relative change in the ratio is assigned as the 
systematic error. The error due to vertexing is obtained 
by taking the difference in the ratio with and without the 
vertex requirement in the D s candidate selection. A sum- 
mary of the systematic uncertainties in signal efficiency 
is given in Table III. Using the measured signal efficiency 



TABLE III: Summary of systematic uncertainties for signal 
efficiencies. 



Systematics 


D-Df (%) D* a -Df (%) D* 3 - 


D* s + (%) 


Tracking eff. 


4.3 


4.3 


4.3 


K s eff. 


2.7 


2.7 


2.7 


Kaon PID 


9.2 


9.2 


9.2 


Photon eff. 




1.8 


3.6 


B counting 


1.1 


1.1 


1.1 


MC statistics 


2.0 


3.5 


4.5 


D ( * } b.f. 


26.0 


26.0 


26.0 


Selection 


5.4 


5.4 


6.0 


Total 


28.7 


28.8 


29.3 



(J2i € iBi)> 211 fb~ on-resonance data corresponding to 
N B g = (231.8 ± 2.6) xlO 6 , the background estimation 
along with the uncertainties and the observed candidate 
events in the signal region N can d, we determine the 90% 
confidence-level (C.L.) upper limit using the procedure 
given in [16]. The systematic uncertainties are included 
following the prescription in Ref. [17]. In all branching 
fraction calculations we assume equal production of B°B° 
and B+B- pairs at the Y(AS). 

The search for B a -> D~ D+ , B° -> D*~D+ and 
B" -> D*~D*+ decays yields the 90% C.L. upper lim- 
its (Table II): 

B(B° -> D~Df) < 1.0 x 10~ 4 , 
B(B° -» D*~D+) < 1.3 x 10~ 4 , 
B(B° -» D*-D*+) < 2.4 x 10~ 4 . 

In conclusion, we have performed a measurement of the 
decay rates for B° -> DjD+, B a -> D*~D+ and 
B° — > D*~ D* s + processes with a sensitivity needed to 
test the SM prediction [18] . Our upper limits disfavor the 
branching fraction predictions in Ref. [3] for all three B° 



decays and accommodate the predictions of the pQCD 
calculation [2] for all three B° decay modes. The possi- 
ble existence of a significant W-cxchange component in 
B° — > D~D + [19] decays is not confirmed in this analy- 
sis. 
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